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The velocity potential, lift force, mioB^nt, and propulsive 
force on a two-dimensional airfoil in a stream of periodically 
varyinff angle of attack have "been derived on the 'basis of non- 
statlonary incompressible potential-flow theory which Includes 
the effect of the continuous sheet of vortices shed from the 
trailing edge . Application of these restilts was made in an, 
analysis of the variation with frequency of the propulsive force 
on an airfoil in an oscillating stream and in an analysis of the 
problem of forced vibrations of an airfoil in an oscillating 
stream with consideration of the stiffness of the airfoil and 
the position of its torsion axis ^ It was shown that -vdien the 
torsion axis of the airfoil is ahead, of the quarter- choi*d point 
the amplitude of vibrations is generally not large, but when the 
torsion axis is behind the quarter-chord point certain conditions 
exist -under which dangerous amplitudes of vibration m^ occur. 

The nonuniform response which was found for a f reply hinged 
airfoil restricts -the use of such a device as a flow-measTiring 
instrument to the measurement of only very low-frequency angular 
variations in an oscillating stream. 

It Is expected that the results of the theoretical trea-tment 
of the propulsive force will be useful in considerations of 
counterrotating-propeller efficiencies and -fchat the analysis of 
the problem of forced vibrations will be usefiil in design con- 
siderations of lifting surfaces operating in oscillating streams; 
for example, -wind-t-unnel fan blades behind a set of prerqtation 
vanes, 'or tail surfaces In fluctuating wakes. 


INTRODUCTION 


The phenomenon of an airfoil in an pscillatlng stream (that 
is, a stream of which the angle of attack varies periodically) is 
encountered in many phases of aeronautics. For example, the 
effect of a set of prerotation vanes upon a wind-tunnel fan blade 
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is to produce periodic distur'bances throu^ vhlch the fan hlades 
pasa. Further, a conalderatlon of counterrotating propellers 
shows that the rearward blades operating in the helical waJce 
produced by the fonrarcl. blades are in a stream of varying angle 
of attach. As another example, the horizontal tail of an air- 
plane may be subjected to forces induced by fluctuations in the 
angle of the wing wahe . - 

The theory of nonstationary motion around airfoils with 
consideration of partial motions of the fluid has been developed 
by Garrick (reference 1), Kttssner (reference 2), and others. The 
present alternative treatment leads to a value for the lift force 
on an airfoil in an oscillating stream which is in agreement 
with that obtained in references 1 and 2. The treatment in 
reference 2 leads to general results for the propifLsive force. 

The present paper, employing somewhat different derivations, 
gives explicit results for the propulsion as well as for the 
lift on an airfoil in an oscillating stream. 

Two problems which arise in cases of oscillating flows 
are (l) the production of vibrations and ( 2 ) the so-called 
"Katzmayr effect" (reference 3 ) or existence of a propulsive 
force. With regard 'to problem (1) the object of the present 
paper is- to examine theoretically the dynamics of an airfoil in 
an oscillating ptroam and to deteymlne unde_r what conditions 
dangerous amplitudes of vibration may occur. A special case 
for which the torsional stiffness is zero is treated with a 
view to the possibility of using a small freely hinged airfoil 
as a device for measuring the angular amplitude of .an oscillating 
stream. Previous work on this problem (reference it) has been 
done for' the case in which the stiffness of the airfoil was 
expected to give large vibrations . With regard to problem (2), 
which is of importance in considerations of counterrotating- 
propeller efficiencies, a theoretical investigation is made of 
the horizontal forces experienced by an airfoil in an oscillating 
stream. 

The theoretical development is divided into three parts: 

( 1 ) derivation of the lift force and moment acting on an airfoil 
in an oscillating stream, ( 2 ) derivation of the propvilslve force, 
and ( 3 ) derivation and solution of the equation of motion of an 
airfoil executing torsional vibrations in an oscillating stream. 
The theoretical methods used iri the derivation of the lift forces 
and moments consist in an extension of the methods of Theodorsen ■ 
(reference 5). For the derivation of the propulsive force 
application is made of the method outlined by von and 

Burgers (reference 6 pp. 52 and 30^)* The following visual assump- 
tions are made throu^out; (a) incompressible potential flow. 
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(■fa) two-dimensional flat-plate airfoil, (c)- ecmall oscillations, 
and (d) plane wake extending frcaa trailing edge to infinity. 


SYMBOLS 


a 

a 

P 

X 

t 

V 

V 
tu 
k 

C(k) 

F,G 

p 

P 

■p 



R. 


lialf chord of airfoil 

x-coordinate of torsion axis of airfoil 

angle of attack of - airfoil ' measured clpckwisevfrom 
horizontal 

angle of alrstream finolSi mean direction- measured 
positive counterclockwise 

horizontal coordinate j nondimens lonal with respect to h 
time 

, stream velocity 
frequency 
circular frequency 
reduced frequeijcy 

Theodorsen’s C-fwnction fi’om reference 5 (F + iO) 
real' and lmaglnary; partB of C— function 
local static ^pressure 

air density _ 

perpendicular force . 

pitching moment about x = a measured positive 
counterclockwise 

moment of inertia, about x = a 

torsional stiffness of wing 

propulsive, .force 
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«P 

'Pr 

U 


noncirculatory velocity potential 
circulatpry velocity potential 
strength of -walce discontinuity 
phase angle 


Jq( k) Bessel fxaictlons of the first kind and zero and 

first order 

To(h),Y2^(k) Bessel fnnctlon of- the second kind and zero and 
first order 


7 

m- 

r„ 


strength of vorticity distrihution on the airfoil 
mass per imit length of .^ing 
radius of gyration divided By h 

ratio of mass of cylinder of air of diameter equal 



to chord of ving to mass of wing 
constant 


SE^ 

m 


LIFT FORCE AND PITCHPIG MOMENT 


In accordance with Theodor sen (reference 5)7 the forces due to 
the noncirculatory flow and. to the effect of the wake are treated 
separately . 

Noncirculatory force and moment .- Consider an airfoil of 
chord 2 b at zero angle of attack with respect to the average 
direction of a sinusoidal stream traveling to the right (fig. l) • 

If the amplitude Pq angle-of-attack change in the stream is 
small then ilie horizontal and vertical components of the velocity, 
respectively, are given By 


and 


■V[j = V cos Pv T • 


vp = V sin p Rf vP 
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whelre, with the assumption of sinusoidal oscillations. 


3 = 


Thus, the airfoil may "be considered as heing in a uniform, 
horizontal stream of velocity . v plus a vertical sinusoidal 
gust of the form 


e(x) = vPoei(^^“^> 


The velocity potential cp for such a normal-velocity distribution 
is (appendix A) 


_ _ (/i ^ Ql(<ot-^a:) 

Ic 


where Jq(u) is a Bessel function of the first kind and zero 
order . 

Use of the equation of motion for nonstationary flow gives 


qIux j^(u) du 


( 1 ) 


3w 

St “ 



where 

w fluid velocity 

p local static pressure 

p air density 
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and the substitution v 
at the point z as 



gives the pressure difference Ap 


(i I 



(2) 


Integration of this local pressure difference over the length of 

the airfoil gives as the- total force P (see, appendix B, equation (B8)) 


P = 2rtioBv%Qj^(k)e^^t (3) 

where J]_(k) is a Bessel function of the first kind and^ first 
order . 

The nonoiroulatory moment about x » a (flg» 2) is obtained 
from the integral 


M. 


r ^ 

b^ / Ap(x - a) 
^- 1 , 


dx 


which yields (appendix B, equation ,(Bl6)) 


\ =* -itpb^v^PQe^'^^laJ-j^(k) + Jo(k)J 


(M 


Circulatojry force and moment .- The 'velocity potential of an 
element of vorticity -AT at a position xq in the wake and its 
mate AP distributed over the airfoil is (reference 5) 
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The element -AP moves to idle ri^t relative to the airfoil with 
a velocity -v. Thus, 


^3CTQ 

St 



SXq 


Substituting this e3T)ression and 


3cp 




into eq.uation (2) and 


Integrating the effect of the entire wake on the airfoil yields the 
force 


P = -pvh / TI dXf. (6) 

- 1 


where U dxQ ie the element of vorticity AF at the point Xq. 

The Kutta condition requires that at the trailing edge of 
the plate the induced velocity equeLlszero; therefore, at x = 1 



= 0 


where 


CPp » -b ^ Cpx3^ dXQ 

Introducing the potential cp from equation (l) resilts in 
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Let 


9 = . !!0 r ei” Jo(u) du (8) 

k , -^0 

vhere u is an integrating varla'ble. Equation ( 8 ) 1:60011108 
(appendix C) 


Q = (9) 

Conibining equations (6), (7)^ and (8) and assuming the wake to he 
of the form 

XT = XJ^ei(a>t--laco) , . 

gives 'for the clarculatory force 


P 


-2rtpvhQ 



(10) 


Similarly, the circulatory moment which is obtained from 


is (reference 5) 


t3 -2rtpvh Q 


Ma = h2 J' Ap(x - a) dx 


/ 3Cq -ikxo 

1 Ya + i Wl ^ 

2' \ 2jji^ Xq + 1 -ikxo 


dXQ 


- 1 


ikxo 

e dxo 


( 11 ) 
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The integral expressi-.n In equations (lO) and (ll) is Theodorsen's C- 
f unction (reference “5) . Thus, 


P = 2rtphv2[3QG(k) ■ Jo(k) - iJ^Ck)] (12) 


and 


Ma = 2irpv2b%Q | - + 0C(k) 


jQ(k) - iJ3^(k) 


.lost 


(13) 


Adding equations (3) and (12) gives for the total force 


r 

P » 2irphv%Q <jC(k) |jQ(k) * lJj_(k)J + iJ^(k) 


,io)t 




This expression agrees with that given "by Garrick (reference 1) and 
Kussner (reference 2) in 'which soffis'what different methods of 
derivation are used from those used in the. present paper. 

Adding equations (4) and (13 ) gives, for the total moment 
about X = a. 


M 


= -irph^v^pQ j2iaJ3^(k) + Jo(k)J 


,io3t 


2irpb^v^3( 




C(k) 


Jq(1c) - iJi(k) 


,iCDt 


(15) 


Examination of equations (lif) and ( 15 ) leads to 


Ma = -■b(a + i 




This equation means that the center of pressure is at the quarter- 
chord point of the airfoil.. 



10 


NACA TN No. 1372 


PROPULSIVE POECE 

According to von Karro^-aiid Btirgers (reference 6, p. 52), if 
the strength of the vortioity distrihution at the leading edge of 
the plate is of ‘the form 




(16) 


then the suction or propulsive foroe acting on the airfoil s . has the 
value 


. ' 8=‘jfpC^ : - ( 17 ) 

vhere C is a constant. The vortex strength is given hy the pvan 
of the tangential velocities on the two sides,.; ' Thus . 



and, therefore. 


C =! 


-<_x 


-I 




(18) 


Carrying out the indicated calculations in appendix D leads to 


■I 


2 0 


ICJit 


+ IJ3J + (Jo - iJi) 


Jl - •!^0 - ^1) 




Ji + V* i(Jo -"^0) 





MCA TN No. 1372 


11 


vhere-^ means JqC^) and so forth and YQ(k) and Y^Ck) are 
Bessel functions of the second Mnd and of the zero and first order. 
The real part of this' eq.uation is . 


C <=! /h •vBo(|/(X + Y) cos- mt- - .^(X - Y) sin cut 




vhere 


X V (A. 



(jo - ’fl)® + (»1 + Yq)‘ 



(Jq - - (Jj * 

(jo - ('ll + 


and finally 


B = (X f Y cob 2o)t - Z sin 2cnt) 

V -i , ■ 

where 


( 19 ) 
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DYNAMICS OF AN AIRFOIL IN AN OSCILLATING STREAM 


The problem treated herein is that of an airfoil executing 
torsional vibrations In an oscillating stream. Because of assump- 
tions (c) and (d) of the Introduction - namely, that the oscillations 
are small and that the wake is plane - the aerodynamic moments due 
to partial motions of the fluid and those due ^ motions of the 
airfoil' can be treated independently . 

The pitching moment acting on an airfoil undergoing angular 
oscillations in a uniform stream is (reference 5) 



Thus, the total aerodynamic moment acting on 'an oscillating 
airfoil in an oscillating stream is the svim of eq.uations (I5) 
and (20). The equation of motion is obtained by expressing the 
equilibrium of the aerodynami.c moment, the moment of inertia, 
and the mechanical restoring moment. Thus, if structural damping 
is neglected, the equation of motion is 


Aerodynamic moment =» IgC^ + 

where 

I^ moment of inertia about x => a 
Ea ^ torsional stiffness . - 

Combining equations (15), (20), and (21) gives 


[ja^ + + a2)J a + 1.2 (1 + ^C(k) 


a 


_ 2kv^ 
b^ 


(a + |)3o ^(^^)(jo - i^i) + 


,i(ot 


( 21 ) 


( 22 ) 



•where 


m mass per unit length of wing 
r„ radius of gyration divided by b 



ib2 , 


mo 


ratio of anass of cylinder of air of diameter equal to chord of wing to mass of •wing 


{¥) 


The differential equation ( 22 ) Is an equation for forced vibrations^ the solution ■to •^rtiich is' 
fotmd by letting 

a s* ' ; (23) 

•v*ere 


oq amplitude of angular oscillations of the wing 


ijf phase' factor 

The ratio of the amplitude of airfoil oecillations to the ampli^fcude of stream oscillations 
•will be, called the response of the airfoil. If the right-hanji side of equation (23) is 
substituted into equation (22), after some-what lengthy but strai^tforv/ard calculations •the 
square of the response of the airfoil la found to be 



where F(lc) + lG(]s:) 


C(k)j Theodorsen ' s C -function. 


( 24 ) 


H 

00 
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DISCUSSION OF liSSULTS 


P ropulsive for ce.- It has "been shown (see eq.i:iatlon (I9)) that 
the propulsive force acting on a fixed airfoil at zero .englp of 
attack in an oscillating stream is 


• s e itp'bv2pQ^(X'+ Y cos 2ti>t - Z Din 2ajt) 


Since Srtphv^Pg is the .stationary vali’e of the lift on on airfoil 
at angle of attack Pq, let 


Lq «= 2fl:pbv^ 


Then 


B s! 


loPo 


(X t Y cos Opt - Z sin apt) 


(25) 


and the average value of s is 




' 0^0 


In figures 3 and :4 ciATves are presented that show the variation 
with wave lengths of stream oscillations of the coefficients X, Y, 
and Z apnearjjig in the aquation for the propulsive force. For very 
lew frequencies ” that is, for long wave lengths of stream oscilla- 
tions - X and y approach 'the value 1, and Z hecomes zero. Thus, 
as k—?0, equation (25) becomes 


• s w 


O' o 
2 


(1 + COB 2 


(26) 


This resxolt is exactly that which is "to be expected from quasi- 
stationary considerations in which the lift is assumed to be 
instantaneously that value prescribed by the geometrical angle of 
attack^ that is, the shed vortlcity produced by variations in 
angle of attack is assumed to appear instantaneously at a point 
infinitely distant from the airfoil. Thus (see fig. 5) the 
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perpendicular force L on an airfoil at angle of attack ^ with 
respect to the stream is 


L = 2rtpv^‘b3 


and hecauae from well-known considerations of two-dimensional 
flows there can he no induced drag, the propulsive force must he 


a = L3 = 2itphv^3® 


and if 3 = 3 q ■ 


a s 23rphv^3o^coB^ajt 

= Lq3o(^| + J 


which is the same as equation (^). 

The reason for plotting the ^coefficients Y and Z against k 
as well as against where is the wave length of stream 


oscillations, is that in the. neighborhood of ~ 

■ iC 


0 the values of 


y and Z fluctuate infinitely many times. This behavior is 
brought' out in fi^re 4 where it may be seen that between . 


k ;= 1 


I- 


and k =03 (i = o) the curves oscillate about 


zero with decreasing amplitvide. 


Forced vibrations .- An aneilysis is made of the response (see 
equation .( 2ij-)) of ah airfoil in an oscillating stream with particular 
emphasis on the parameters a and Sg^, these parameters having 
qualitative as well as quantitative effects on the values of the 
response. .The stiffness in torsion Rg^ is related to the natural 
frequency V ' for zero stream velocity by the equation 



i6 


MCA TN No, 1372 



where X_ 

cl 

airfoil . 


+ jrph^^ + is the 

From equation (27), 


total moment of inertia of the 


B 


a 


2Kmv2 




+ 




where k' = 

V 


In figure 6 are shown fospohee curres for on airfoil hinged at 
the leading edge (a = -l.O) and haring various natural frequencies 
of oscillation. The following values of the parameters k and r^ 
were chosen as "being within the practical range of application: 

K = 0.0653; ^ (flat- plate mass distrllJUtion) . The response 

for any natural frequency is never very large iat most a"bout — » 2.5]* 

\ 3o / 

and the response decreases with Increasing stiffness. The 
value of the stream frequency at which maximum response occurs is 
seen to correspond more closely to V as the natural frequency 
increases. Even for a freely hinged airfoil (k' ■= 0) ' a sort of . 
resonance frequency exists. (See -fig. J.) Thus the use. of such 
a device for measuring angular variations in an oscillating stream 
would "be valid only in the range of very low reduced frequency (long 
wave lengths)' in which the response approaches unity. 

Somewhat different phenomena occur when the hinge point is 
"behind the quarter- chord point of the airfolD, . In general it may 
"be stated that the response is gi’eater. In particular, q. critical 
stiffness exists "below which the airfoil is in unsta"ble equill"briUDi. 
The condition for divergent motions of the airfoil is that the 
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coefficient of a in equation ( 22 ) "be less than or equal to 0 
for k =s 0 . The criticAL stiffness, therefore, is defined, hy 




2Kntv^ 


= a 




In figure 8 are presented some results for the hinge placed 
near the center of the airfoil. The values of the parameters a, 
K, and ra^ chosen -were; a = -0.1, k = o.06‘?3, and rg^^ = 

The reduced critical frequency foi’ this case is k' = O.385. For 
stiffness values somewhat hi{dier than the critical, the response 
is not unduly large and again, as when the hinge was at the 
leading edge, the maximum response with large values of the 
stiffness occurs at a stream freqtiency close to the frequency v ' , 


CONCLUSIONS 


The velocity potential, lift force, moment, and propulsive 
force on a two-dimensional airfoil in a stream of periodically 
varying angle of attack has been derived on the basis of non- 
stationary incompressible potential-flow theory which includes 
the effect of the continuous sheet of vortices shed from the 
trailing edge. Application of these results was made in an 
analysis of the variation with frequency of the propulsive forces 
on an airfoil in an oscillating stream and in an analysis of 
the problem of forced vibration of an airfoil in an oscillating 
stream with consideration of the stiffness of the airfoil and 
the position of its torsion axis. Ilie following conclusions 
were indicated: 

1 . The value of the propulsive force acting on an airfoil 

in an oscillating stream is sufficiently large to be of practical 
importance . 

2 . The amplitude of vibration of an airfoil in an oscillating 
stream is critically dependent on the stiffness of the airfoil 
and the position of its torsion axis. In general, amplitudes of 
vibration are smaller when the torsion axis is ahead of the 
quarter-chord point and larger when the torsion axis is behind 
the quarder-chord point. Because of the nonun.lform response of 
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a freely hlnfted airfoil the use of such a, device for the measurement 
of angular variations in an oscillating stream would he restricted 
to the range of very low frequency in which the response approaches 
unity. 


Langley fiemorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., Apx'il 28, 194-7 
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APPENDIX A 

NONCIRCnLATOEr TELOCITY POOEHTIAL FOR M 
AIRFOIL IN A SINUSOIDAL GUST 


The prohlem of finding the velocity potential for an airfoil 
having a certain normal velocity dis trihut ion is solved hy the 
method Indicated hy von Kairmsri and Burgers (reference 6, p. 44) . 

Represent iiie wing hy a circle (fig* 9). Place a source of 
strength 2e at the point y^) on the circle and a sink of 

strength -26 at . The velocity potential of this 

source- sink pair in the notation of the present paper is given hy 
(reference 5) 


Ap = 


6 

2 It 


(x - x^)^ + (y - yi)2 

log - — I i 

{x “ Xj}® +-(y + y^)^ 


The transformation of the circle to its diameter is 


y = ^1 - x^, X = X 
For the airfoil in a sinusoidal gust 


6 = 


Thus 


Ap = e^'^e'i^l 

2rt 


+ (y - yi)^ 


(x - xi)2 + (y + yi)^ 



20 


WACA TK No. 1372 



Integration ty parts leads to 
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But, from relation (U) on pa^o of reference 7, 



vhere Jo(k) is a Bessel function of the first kind and zero order. 
Therefore, 


^ = irtjQ(k) - ixf (A2) 

Equation (A2) is a non-homogeneoup differential equation of the first 
order. The homogeneous part 


8f 


-ixf 


has as solution 


f = 


ce 


-ikx 


where c is an arhitrary constant with respect to k. The 
particular solution is obtained by the meth^ of variation of 
parameters ( reference 8, p . Ilk) , Let 


f e g(k, x)e~^^ 


(A3) 


-ikx 


Then 


M ^ e-i^ - ixge 
Sk 5k ^ 
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Combining this expression with eq.uatlon (A2) gives 

|8 = 

integration of ■which leads to 

g = irt r e^^Jrt(u) du (A4) 

Comhlning equations (A3) and (a 4) gives 

s^^Jq(u) du 

which, when substituted into equation (Al), gives for the non- 
circulatory velocity potential 


f(k; , x) » Irto"^^ 


r 
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APPENDIX B 

LUT AND ^^0^^ENT FOR WONCIRCULATORy FLOW 
Lift 


n?he total force on tBe airfoil is given "by 




Ap d;!C 


vhere, from equation (2) 


Ap = -2p |Z 1$ + ^ 

vb Sx at 


But 


Therefore 


'Px-l = 'Px*-! “ ^ 


P = -2ph 


f ^ ^ 

“^-1 at 


dx 


From eq.mtion (l). 


I = . /l . f e‘“j„{u) 


du 


(Bl) 


Thus 


P 


i2ph^vBcP 


aiot 


k 


\/l - dx J' 6^^Jq(u) du (B2) 
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Let 


fl(k) « y' /l - dx y* e^^JgCu) du (B3) 


Interchanging the order of integration gives 


fj^(fc) a J' Jq<u) du ly 


2 l(u-»k)x 


X e 


dx 


Prom reference 7> page ^8, 


jv(.) = yy.i\ r 

■■ C'’ ♦ £)<I7 


Therefore 


(Bl;) 


(B5) 


j (u - k) « ■ ^} . r ^ (1 - x2)l dx 

. a- i2)i ai • 

TT 


Suhstitutlng this erpreseion into equation (BU) gives 


'fl(k) a..« Jo(u) - k> 
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Let 


and 


Then 


u - k = -V 


du B -dw 


fide) = ,t J Jo{k - w) s 

k 


dw 


= -ir r : J,(-v)'Jw(k - w) ^ . 

v/q .4. U y 




J'oO' - ■') ? 


But from reference 7, page 380 , 'for u > 0 and v.> -1 


yt) J^(Z . t) 


(BS) 


Therefore 


fl(k) B rtj^(k) 


(B7) 
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and equation (B 2 ) 'becomes 




(B8) 


Moment 

The nonolrciilatory moment about the point x * a (fig. 2 ) 
is o'btained from the integral 


Ma " 


J' Ap(x - a) dx 
which, combined with eqtiation ( 2 ), gives 


Mg, s 2 pbv J' cp dx + 2 pb^a J' ^ dx - 2pb2 J' 




By use of equations (l) and (Bl), this equation becomes 


2pb^v^pf. pi — ^ fic 

= - t/l - x 2 dx J^(vl) du 




,iux 


. gipgWgQ r ^ 


k 


/ </“ 

^-1 


^ e”^^ dx j e^"^ Jo(t^) (B 9 ) 


/ 


iux 


2ipb-^ojvPQ 


.lOJt 




r 


^2 0-ikx / giux J (t;i) a.u 
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The first two integrals on the rl^t side ^of this expression are 
already known (see equations (B3) snd (B?))* In order to obtain 
the third integral, let 


fo(k) o r X /l - dx f 6^^ (bio) 

sZi ^0 


Interchange the order of integration. Then 





Jo(n) du 



gi(u-k)x 


dx 


(Bll) 


But 



^i(u-k)x ^ 


i'(l - x2)3/2 gi(u-k)x 


1 

-1 


(B12) 


+ r (1 - to 

« j;) f (1 - x2)3/2 ^±(n-Ts:)x 

3 - • -1 


and, from equation (B5) 



, ^2)3/2 3i(u-k)x ^ ^ :■ 

(u - -k)2 


(B13) 
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Combine equations (B12) and (B13) end substitute into '^qiiation (Bll) . 
Then 


f^Ck) » iJt r Jq(u) J2(u - k) — ^ 


Let 


and 


Then 


u - k «s -V 


du = -dv 


fgCk) = y' "^0^^ - w) ^ <= -iJt y JgC-w) J^Ck - w) ~ 


» ,-i« JgCv) Jo(k - w) ^ 


which hecomes (see equation (b6)) 


fp(fc) - - |2 J,(k) 


(Bl4) 


The equation for the moment about x *» a can now be written as 




2Ttpb^^Po . 23Tipb^ao>v|3o , ^ 

- -J: £ Jl(k) e^^ r ^ J (k)e^^t 


k 

rtpb3MVpQ 


(BI5) 


J2(k) e^^^ 
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The recurrence formula from reference 7, page 17, 


’^n-lC^^ ^ ^n(^) 

gives for n « 1 

J2(k) » -JoC^^) + I Jl(^) 


Substituting this expression into equation (BI 5 ) and making use of 
the definition k = ^ yields 


= -2rtiph2aT2pQ jf^(k) e^^.-;-«pB^“13oJo(k) 



APPENDIX C 


EVALUATION . OP 0 . (EQUATION • ( 8) ) 


Equation (8) is 


Q = - 2° r Jq(u) du 

k 'Jq V 

^0 ia>t i(u-k) _ . V - 

s» - e / ® 

1: Jq 


But 




cos (u - k) Jq(u) du 


+ i 



Bin (u - fc) Jq(u) du 


From reference J, pages 380 and 381, 


COB (k - u) Jq(u) du = 


J sin (k - u) Jq(u) du = kJ 2 (k) 
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Therefore 


y* ji(ii-k) 


Jq(u) du = fcTo(lc) - IkJi(k) 


= kjjo(k) - Uj^Ck)'] 


and 
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APPENDIX D 

CALCULATION OF THE PROPULSIVE FORCl 
From equation (l) 




^^0 Ql(<iit+k) 
k 




e"^“ J (u) du 




The Integral on the right aide of thia equation ia the complex 
conjugate of that in appendix C. Thua 


Jg(u) du = k(jg - ijj^) 
where Jq meana Jg(k), and so forth. Therefore 




Now, for the velocity dvi© to the wake. 


.icut 


(Dl) 
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tut (see equations (7) and (8)) 



From reference 7, page l8o, 


sin cosh t - i n cosh nt dt 


and 


Y (k) = - ^ / cos (k cosh t 


^ cosh t - i n it^ 


cosh nt dt 
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vhich lead to, vhen cosh t Is replaced ‘by 3 Cq, 


Jo(k) 


2 


sin k xq dX (3 




1 


Jl(k) 



XqOOB k'XQ dXQ 



Yo(&) « 



cos k Xq dXQ 

i/xq2 - 1 


Yi(k) 



XqSIh k Xq dXQ 
V^Xq2 - 1 


Equation (D2) hecomes now 



„ (Ji - * W 

C*^i ^o) 


and after substitution of the value of Q from equation (C2) 
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= - f ^ (Jo “ iJl) 

‘px / 3>=-l \j/l + ^x=-l i/2 


(^Jl - Yq) - i(Jo + yJ 


.iojt 


}_(Jl + Yo) + i(Jo - Yi)j 
^ ' (^3) 

Coiabining equations (Dl) and '(D 3 ) ' 5 ^th“(see equation (I8)) 


C = 






gives 


-1/1 (-3b + iJl) + (.Jo - IJl) ~ 

> ^ /' T- J. V- \ a. ^ Z' _ V, \ 


r .. 


(^Jl + Yo) + l(Jo - Yi) 


- tPq(cos cut + i sin <x)t)'<(Jo 




1 + 




K . 


gJl(jQJl YqYi) 

(Jl + Yof + (Jo - Yi)‘ 


ir iJi 


(Jl + Yo )2 + (Jo - Yx)' 


Jl^ - Jq^ + Yl^ - Ypg 
(Jl + Yo )2 + (Jo - Yi )2 


_i 2 Jo(JqJi YqYi) [ 

(Jl + Yo) 2 + (Jo - Y]^)2 I 

the real part of yhich Is 


C =t /b vp 


0 


l/^X + Y) cos cot - \/(X - Y) sin cot : 


(335) 
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where 


(Jo - Tl)® ♦ (Jl + '^of 


. — i_ 

and 

Y „ i(J 0 - - (Jl ^0)^1 

f(Jo - 

= /J.f (Jq - - (Ji ♦ ^0;^ 

vtk/ r 1 2 

. L'!>'^o - * (Ji + ^o;®J ■ 

and use has "been made of the formula (reference 7) P» 77 ) 


Jo^l - Vo = - 4 


Suhstltuting eq.uation (D 5 ) into (see equation (l?)) 


s =, JtpC^ 
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gives 


s = . itpTsv^Pq® 


(X + Y) cos^t + (X 


-2 



sin mt cos cot 




= itpTsv'^Po'^ fX + Y cos .2 03t - ^ 


Letting Z = /x^ - Y^ . gives finally 


s 


irpliv^o^. (X + Y cos 2 ost - 


- Y). sin^nt 


- Y^ sin 2 dit 


Z sin 2 cot) 
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